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ABSTRACT 

We investigate peculiar velocities predicted for clusters in Lambda cold dark mat- 
ter (ACDM) models assuming that the initial density fluctuation field is Gaussian. To 
study the non-linear regime, we use N-body simulations. We investigate the rms veloc- 
ity and the probability distribution function of cluster peculiar velocities for different 
cluster masses. To identify clusters in the simulation we use two methods: the stan- 
dard friends-of-friends (FOF) method and the method, where the clusters are defined 
as maxima of a smoothed density field (DM AX). The density field is smoothed with a 
top-hat window, using the smoothing radii R s = 1.5/i -1 Mpc and R s = l.O/i -1 Mpc. 
The peculiar velocity of the DMAX clusters is defined to be the mean peculiar veloc- 
ity of matter within a sphere of the radius R s . We find that the rms velocity of the 
FOF clusters decreases as the cluster mass increases. The rms velocity of the DMAX 
clusters is almost independent of the cluster mass and is well approximated by the 
linear rms peculiar velocity smoothed at the radius R = R s . The velocity distribution 
function of the DMAX clusters is similar to a Gaussian. 

Key words: galaxies: clusters: general - cosmology: theory - dark matter - large- 
scale structure of Universe. 



1 INTRODUCTION 



One of the interesting problems in cosmology is the evolution 
of the large-scale peculiar velocity field in the Universe. The 
evolution of the large-scale velocity field has been studied in 
a number of papers (see, e.g., Davis et al. 1985; Kofman et 
al. 1994; Brainerd & Villumsen 1994; Gramann et al. 1995; 
Jenkins et al. 1998; Juszkiewicz, Springel & Durrer 1999; 
Colin, Klypin & Kravtsov 2000; Gramann & Suhhonenko 
2002; Ciecielag et al. 2003). The linear rms peculiar velocity 
on a given scale R can be expressed as 



a v (R) = H f 



(1) 



where Hq is the Hubble constant and / is the dimension- 
less growth rate, which is related to the cosmological matter 
density parameter, f2 m , and the cosmological constant, Qa, 
by 

/(n m ,n A )~n^ 8 + ^.(i + ^) (2) 

(Lahav et al. 1991). The spectral moments (Jj(R) are defined 
for any integer j by 



o-](R) 



1 

2^2 



P(k)W 2 (kR)k 



2j+2 



dk, 



(3) 



where P(k) is the power spectrum of density fluctuations 
and W(kR) is a Fourier transform of the smoothing window. 



For the top-hat window in real space that we will use in this 
paper, IF (a?) = (3/a; 3 ) [sin(a;) — xcos(x)]. Notice that the 
predicted rms peculiar velocity depends both on cosmology 
and on the shape of the power spectrum. 

Jenkins et al. (1998) used a N-body simulation to inves- 
tigate the rms peculiar velocity of the dark matter smoothed 
on the scales R = 10 — 80/i -1 Mpc. They compared the re- 
sults with the linear approximation (1) and found that on 
scales above 20/i _1 Mpc, the linear theory prediction agrees 
very well with the simulation. Kofman et al. (1994) and 
Ciecielag et al. (2003) studied the evolution of the pecu- 
liar velocity one-point distribution function, beginning with 
Gaussian initial fluctuations. They showed that on mildly 
non-linear scales (4 — 10/i _1 Mpc, Gaussian smoothing) the 
distribution of the Cartesian components of the peculiar ve- 
locity field is well approximated by a Gaussian. On smaller 
scales the velocity distribution of dark matter becomes non- 
Gaussian due to the motions of matter within dense systems 
(see e.g. Sheth & Diaferio (2001) for the discussion of the 
velocity distribution of dark matter particles in different cos- 
mological models). 

Important information for the large-scale velocity field 
is provided by peculiar velocities of galaxy clusters. The 
evolution of peculiar velocities of clusters in different cos- 
mological models has been examined in several papers (e.g. 
Bahcall, Gramann & Cen 1994; Croft & Efstathiou 1994; 
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Moscardini et al. 1996; Suhhonenko & Gramann 1999; Col- 
berg et al. 2000; Sheth & Diaferio 2001; Suhhonenko & Gra- 
mann 2003; Hamana et al. 2003). Clusters represent high- 
density maxima in the dark matter density field. Bardeen 
et al. (1986) studied the peculiar velocity distribution for 
peaks of a Gaussian density field. They showed that the ve- 
locity distribution function for peaks is also Gaussian, with 
the rms peculiar velocity 

Op ( R) = a v (R) y/l- a*/v*a*_ ~. (4) 

Peaks have lower peculiar velocities than field points. The 
reason for this difference is the fact that the velocity field is 
correlated with the density gradient field. However, at radii 
R ~ 1 — 5ft -1 Mpc, the difference between a v (R) and u p (R) 
is small in standard cosmological models. For example, in the 
flat Q, m = 0.3 ACDM model, it is about 3 per cent at the 
radius R — 1.5ft _1 Mpc. Gramann et al. (1995) studied the 
large-scale velocity and potential maps of clusters and of the 
matter using N-body simulations. They found that the large- 
scale velocity and potential fields are recovered remarkably 
well by using the velocities of clusters, 

In this paper we study peculiar velocities of clusters in 
a ACDM model. We use the N-body simulation carried out 
by the Virgo Consortium for the flat ACDM model with 
the density parameter fi m = 0.3 (Jenkins et al. 1998). To 
identify clusters in the simulation we use two methods: (1) 
the standard friends-of-friends (FOF) method and (2) the 
method, where the clusters are defined as the maxima of the 
smoothed density field (DM AX). We use a top-hat window 
with the smoothing radii R s = 1.5ft _1 Mpc and R s = lft -1 
Mpc. The velocity and the mass of FOF clusters are defined 
to be as mean velocity and the total mass of all the particles 
in the cluster. To determine the mass and velocity of DMAX 
clusters, we use the same smoothing as for the density field. 
We study the rms velocity of clusters for different cluster 
masses. The rms peculiar velocity of clusters was studied 
also by Suhhonenko & Gramann (2003). They used a some- 
what different DMAX method. The densities and peculiar 
velocities of clusters were determined using the cloud-in-cell 
(CIC) scheme. In this paper we use a top-hat window, which 
allows us to simply calculate cluster masses. 

We compare the rms velocity of clusters in the simula- 
tion with the linear approximation (4) . What window func- 
tion, W(kR), and smoothing length, R, for the spectral mo- 
ments in the equation (4) are appropriate for cluster veloci- 
ties? Different studies have used different window functions 
and smoothing lengths. Croft & Efstathiou (1994) and Bah- 
call, Gramann & Cen (1994) used Gaussian smoothing with 
a radius R = 3ft -1 Mpc. Suhhonenko & Gramann (1999) 
adopted top-hat smoothing with a radius R = 1.5ft~ Mpc. 
To compare model predictions with the observed peculiar 
velocities of galaxy clusters, Borgani et al. (2000) used Gaus- 
sian smoothing with a radius R — 1.5ft -1 Mpc. Colberg et al. 
(2000) and Sheth & Diaferio (2001) connected the smooth- 
ing length R in the equation (4) with the cluster mass M. 
They studied the rms peculiar velocity, a p (R), at the linear 
radius R = Rl{M), where the linear radius Rl is defined as 
M = Aiv/SpbRt (here pb is the mean background density). 
Similar approach was used recently by Hamana et al. (2003) . 

The velocity of DMAX clusters is determined as the 
mean velocity of matter within a sphere of radius R s = 
1.5/i _1 Mpc (or R s = lft _1 Mpc). In other words, we use 



top-hat smoothing with a radius R s . For this reason, we 
adopt also the top-hat window function and the smoothing 
radius R = R s in the equation (4). We consider also the rms 
peculiar velocity, a p , at the radius R = Rl(M). 

The sizes of FOF clusters are not fixed. If we use the 
FOF method, the mean size of high-mass clusters is larger 
than the mean size of low-mass clusters. Therefore, it is not 
possible to connect the rms velocities of FOF clusters with 
a single radius R — R s . We can use the approximation 
R = R L (M). Note that for the DMAX clusters, the ratio 
Rl/R A s = PI ' Pb, where p is the mean cluster density within 
a sphere of radius R s - For p/pt w 200, the ratio Rl/ Rs ~ 6. 

We also investigate the velocity distribution function 
of clusters. We examine the distribution of one-dimensional 
peculiar velocities for different cluster masses. Do the clus- 
ters exhibit a Gaussian distribution of peculiar velocities? In 
addition to studying cluster velocities at the present time, 
we also study how cluster velocities evolve from some early 
time to the present. 

This paper is organized as follows. In Section 2 we de- 
scribe the algorithms that have been used to identify clusters 
and to determine their masses and velocities. In Section 3 
we study the rms peculiar velocity of clusters for different 
cluster masses and compare the results with a p (R) at dif- 
ferent radii. The velocity distribution function is analyzed 
in Section 4. In Section 5 we briefly discuss the evolution of 
cluster velocities. A summary and discussion are presented 
in Section 6. 



2 SIMULATION OF CLUSTERS 

We study cluster velocities in a N-body simulation carried 
out by the Virgo consortium for the flat ACDM model with 
a cosmological constant. These simulations are described in 
detail by Jenkins et al. (1998). The simulations were created 
using an adaptive particle-particle/particle-mesh (AP 3 M) 
code as described by Couchman, Thomas & Pearce (1995) 
and Pearce & Couchman (1997). In the Q, m = 0.3 ACDM 
model studied here, the power spectrum of the initial con- 
ditions was chosen to be in the form given by Bond & Efs- 
tathiou (1984), 

P<yV> = [l + {aq+(bq) 3 / 2 + (cq) 2 Yf'^ ^ 

where q = k/T, a = 6.4ft _1 Mpc, b = 3/i _1 Mpc, c = 
1.7ft -1 Mpc, v = 1.13 and T = Q h = 0.21. The normal- 
ization constant, A, was chosen by fixing the value of as 
(the linearly extrapolated mass fluctuation in spheres of ra- 
dius 8/i _1 Mpc) to be 0.9. The initial density fluctuation field 
was assumed to be Gaussian. 

The evolution of particles was followed in the comoving 
box of size L — 239. 5ft ~ Mpc. The number of particles was 
N p = 256 3 . Therefore, the mean particle separation \ p = 
L/N p /3 = 0.9355ft" 1 Mpc and the mass of a particle m p = 
p b \ 3 p = 6.82 x 10 10 ft _1 M Q . In the flat Q m = 0.3 model, the 
dimensionless growth rate / = 0.513 (Peebles 1984). 

We used two different algorithms to identify clusters in 
simulation: the standard friends-of-friends (FOF) algorithm, 
and the algorithm, where clusters are defined as maxima of 
the smoothed density field (DMAX). 
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The friends-of-friends group finder algorithm was ap- 
plied using the program suite developed by the cosmology 
group in the University of Washington. These programs are 
available at http://www-hpcc.astro.washington.edu. To test 
our FOF output data, Suhhoncnko & Gramann (2003) in- 
vestigated the mass function of clusters. The cluster mass 
function in the Virgo simulations has been studied in detail 
by Jenkins et al. (2001). Suhhonenko & Gramann (2003) 
found that the agreement between our results and these ob- 
tained by Jenkins et al. (2001) is very good. 

The FOF cluster finder depends on one parameter 6, 
which defines the linking length as b\ p . The conventional 
choice for this parameter is b = 0.2 (see e.g. Gotz, Huchra 
& Brandenberger 1998; Jenkins et al. 2001). In this paper 
we also define clusters by using the value b = 0.2. We also 
study velocities of the clusters defined by the parameters 
b = 0.15. In the limit of very large numbers of particles per 
object, FOF approximately selects the matter enclosed by 
an isodensity contour at pb/b 3 . 

We studied FOF clusters that contained at least ten 
particles. The three-dimensional peculiar velocity of each 
cluster was defined as 



We studied the rms density contrast and the rms pecu- 
liar velocity on the 256' grid. The rms density contrast on 
the grid was 5.05 and 7.29 for the radii R s — 1.5ft _1 Mpc and 
R s — 1.0ft -1 Mpc, respectively. The rms peculiar velocity, 
o~ v , was determined for the fraction of grid points, F, where 
the number of particles N > 1. If there are no particles in 
the neighbourhood of a grid point, the velocity field is un- 
determined. For R s — 1.5ft -1 Mpc, we found that F = 0.96 
and a v = 473 kms" 1 . For R s = 1.0ft -1 Mpc, F = 0.69 and 
a v = 481 kms -1 , respectively. Here we took into account 
the finite size of the simulation box (see next section). We 
also studied the rms density contrast and the rms peculiar 
velocity on the 350 3 grid and found similar results. 

For each DMAX cluster, we investigated the cluster 
mass, M, and the peculiar velocity, v c i, at the radius R s . 
The mass in the cluster was determined as M — N d m p , 
where Nd is the number of particles in a sphere of radius 
R s around the centre of the cluster. The peculiar velocity of 
each cluster was defined as 



N d 



(9) 



Vd = 



1 



1=1 



(6) 



where Nf is the number of particles in the cluster and Vi is 
the peculiar velocity of the particle i in the cluster. 

We also selected DMAX clusters using the following 
method. 

(1) We calculated the density contrast on a grid. For 
each grid point, the density contrast was determined as 



where Vi is the peculiar velocity of the particle i in the 
DMAX cluster. 

Suhhonenko & Gramann (2003) compared the cluster 
peculiar velocities defined by different methods (i.e. the FOF 
method versus the DMAX method) for massive clusters. By 
using different methods to identify the clusters, we select 
almost the same objects in the simulation. But we assign 
different velocities to the same clusters. 



(7) 



where N is the number of particles in the sphere of radius 
R s around the grid point, and 



L 3 



(8) 



is the mean number of particles in the sphere of radius R s . 

(2) We found the density maxima on the grid. The grid 
point was considered as a density maximum, if its density 
contrast was higher than the density contrast in all 26 neigh- 
bouring grid points. The location of the grid point, where 
the density contrast had a maximum value, was identified 
as the candidate cluster centre. 

(3) The final cluster list was obtained by deleting the 
candidate clusters with lower density contrast in all pairs 
separated by less than the radius R s . 

In this way we define the clusters as maxima of the 
density field smoothed by a top-hat window with a radius 
i? s . The smoothing length sets a lower limit on the size of 
detected structures in the simulations. 

We used the smoothing radii R s = 1.5ft -1 Mpc and 
Rs = 1.0ft -1 Mpc. The mean number of particles in spheres 
of Rs = 1.5/i _1 Mpc and R a = 1.0ft -1 Mpc is N = 17.26 
and N = 5.12, respectively. To select the R s = 1.5ft -1 Mpc 
clusters, we used a 256 3 grid (the cell size I = 0.936ft -1 
Mpc). For Rs = 1.0ft -1 Mpc, we used a 350 3 grid (I = 
0.684ft -1 Mpc). For comparison, for the R s = 1.0ft -1 Mpc 
clusters we used also a 256 3 grid. 



3 THE RMS PECULIAR VELOCITY OF 
CLUSTERS 

To determine the rms peculiar velocities of clusters, we used 
the equation 



2 . 2 
v s + V L 



1 \ 2 2 



(10) 



where the parameter v s describes the dispersion of cluster 
velocities, v c u, derived from the simulations and the param- 
eter vl is the linear contribution from the velocity fluctua- 
tions on scales greater than the size of the simulation box 
L. It is given by 



2 mr 



P(k)dk. 



(11) 



N r j is the number of clusters studied. Using eq. (4), the 
linear rms peculiar velocity of peaks can be written as 



CTp(R) = <Ty{R) - 



jj2 y.2 O-q(R) 



o-KRY 



(12) 



The second term in this expression is not sensitive to the 
amplitude of large-scale fluctuations at wavenumbers k < 
2-k/L. Therefore, the linear rms velocity of peaks can be 
expressed, approximately, as 



Op{R) ~ CTp(R) + v 2 L , 



(13) 
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Table 1. The number of clustc 


rs, N c i , in 


liffcrcnt m 


ass interva 


M 


DMAX 


DMAX 


FOF 


FOF 




Rs = 1.5 


R s = 1.0 


b=0.2 


b=0.15 


10 12 - 5 x 10 12 




70153 


51841 


47056 


5 x 10 12 - 10 13 


10647 


9989 


7326 


6764 


10 13 - 5 x 10 13 


10049 


8866 


6542 


5720 


5 x 10 13 - 10 14 


1261 


1052 


856 


695 


10 14 - 5 x 10 14 


688 


470 


511 


393 


5 x 10 14 - 10 15 


23 


4 


38 


24 



where <j' p (R) is determined by the power spectrum at the 
wavenumbers k > 2n/L and vl is given by eq. (11). For the 
ACDM model studied here, we found that vl = 220 kms -1 . 

As a first step, the one-dimensional distribution of clus- 
ter velocities can be approximated as a Gaussian distribu- 
tion (see next section for the study of the velocity distribu- 
tion of clusters). If the one-dimensional velocities of clusters, 
v xi , follow a Gaussian distribution with a mean v x = and 



a dispersion a , then the sum 



X = -o 



\ ^ 2 

v cli 



1 

i=l 



(14) 



is distributed as a \ 2 distribution with the number of degrees 
of freedom v — 3N c i. In this case, the rms error for the 



Avi 



3Nd 



2 
V.,. 



(15) 



We used eq. (15) to estimate the error bars for the rms 
velocities of clusters. 

The clusters were divided into subgroups according to 
their mass. We studied the rms peculiar velocity of clus- 
ters in subgroups for which the mass was in the ranges 
(10 12 - 5 x 10 12 )/i _1 M Q , (5 x 10 14 - 1O 15 )/i _1 M . For 
the DMAX clusters denned with R 3 = 1.5fe _1 Mpc, we 
considered the clusters with masses M > 5 x 1O 12 /i _1 M 
(N/N > 4.24). Table 1 shows the number of clusters and 
Fig. 1 demonstrates the rms peculiar velocity of clusters in 
different mass intervals. The rms velocities are shown for the 
intervals for which the numbers of clusters are N cl > 10. 

Fig. la and Fig lb demonstrate the rms peculiar veloc- 
ities for the DMAX clusters. Fig. la shows the results for 
the radius R s = 1.5/i _1 Mpc and Fig. lb for the radius 
R 3 = l.Oh- 1 Mpc. We see that the rms velocity of DMAX 
clusters is almost independent of the cluster mass. The rms 
velocity of clusters somewhat increases as the cluster mass 
increases. However, this increase is small (~ 7 per cent). 
For a given scale R s , the rms velocity of clusters is simi- 
lar to the rms velocity for the full grid. The rms velocity 
of R s = 1.5/jr 1 Mpc clusters is 471 kms -1 in the mass 
interval (5 x 10 12 - lO 13 )/i _1 M and 507 kms" 1 in the 
mass interval (5 x 10 14 - 10 15 )/i _1 M Q . The rms velocity 
of R s = l.Oh -1 Mpc clusters is 495 kms -1 in the mass in- 
terval (5 x 10 12 - 10 13 )/i _1 M Q and 510 kms" 1 in the mass 
interval (10 14 - 5 x lO 14 )/i _1 M . The rms velocities for the 
R s = 1.5/i _1 Mpc clusters are smaller than the rms veloci- 
ties for the R s = 1.0h~ Mpc clusters. 

The open circles in Fig lb show the results for the clus- 
ters defined by using the 350 3 grid. For comparison, we plot- 
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Figure 1. The rms peculiar velocities of clusters for different 
cluster masses, (a) The velocities for the DMAX clusters defined 
with R s = 1.5h~ Mpc (open circles). The solid circle represents 
the result obtained by Colberg et al. (2000). The dashed line 
shows the linear rms peculiar velocity of peaks, cr p , for the radius 
R = 1.5ft" 1 Mpc. (b) The velocities for the R s = LOh' 1 Mpc 
DMAX clusters. The circles show the velocities for the clusters 
defined using a 350 3 grid. The dot-dashed line describes the clus- 
ter velocities for a 256 3 grid. For comparison, we show the linear 
rms peculiar velocity of peaks for the radius R = l.Oh^ 1 Mpc 
(dashed line), (c) The velocities of the FOF clusters determined 
by b = 0.2 (circles) and b = 0.15 (triangles). The dotted line in 
each panel shows the linear rms peculiar velocity of peaks for the 
radius _R = R L (M). 
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ted also the rms velocities for the clusters defined by using 
the 256 3 grid. We see that the rms peculiar velocities of 
clusters in different mass intervals are not sensitive to the 
number of cells in the grid that was used to find the cluster 
centres. 

Our results are in good agreement with the results ob- 
tained Colberg et al. (2000). They also studied the rms veloc- 
ity of clusters in the ACDM model, but used a slightly differ- 
ent method to select clusters. High-density regions were lo- 
cated using a FOF method with b = 0.05 and their barycen- 
tres were considered as candidate cluster centers. Any can- 
didate centre for which mass within 1.5ft -1 Mpc exceeded 
the threshold mass Mt was identified as a candidate clus- 
ter. The final cluster list was obtained by deleting the lower 
mass candidate in all pairs separated by less than 1.5ft -1 
Mpc. The peculiar velocity of each cluster was defined to 
be the mean peculiar velocity of all the particles within the 
1.5ft -1 Mpc sphere. This method to determine the cluster 
velocities is similar to the DMAX method, if we use the 
radius R s — 1 .5ft -1 Mpc. 

In Fig. la we show the result obtained by Colberg et 
al. (2000) for the ACDM model. They used the value M t = 
3.5 x lO 14 ft -1 M0. For this value, the number of clusters 
was N c i — 69. They found that the rms cluster velocity 
derived from simulation is v 3 = 439 kms -1 . If we include 
the dispersion v\ , we find that this value of v s corresponds 
to the rms velocity v rm s =491 kms -1 . For comparison, wc 
also studied the R s = 1.5ft -1 Mpc clusters with masses M > 
3.5 x lO 14 /i _1 M0 . We found 72 clusters with the rms velocity 
Vrms = 4821^ kms -1 . This value is in good agreement with 
the rms velocity found by Colberg et al. (200). 

For comparison, we show in Fig. la and Fig. lb the rms 
peculiar velocity of peaks, a p (R), for the radius R — R s . For 
the radii R = 1.5k' 1 Mpc and R = 1.0ft -1 Mpc, the a p = 
500 kms -1 and a p = 510 kms -1 , respectively. We see that 
the rms velocity of DMAX clusters is well approximated by 
the linear rms velocity smoothed at the radius R — R s (i.e. 
at the radius which is used to define the cluster velocity). 
The rms velocity of low-mass clusters is somewhat smaller 
than predicted by the linear theory. 

In Fig. lc we show the rms peculiar velocities of FOF 
clusters determined by 6 = 0.2 and 6 = 0.15. The rms veloc- 
ity of FOF clusters decreases, as the mass of the clusters in- 
creases. This result is in agreement with the results obtained 
by Sheth & Diaferio (2001) and Hamana et al. (2003). They 
studied the rms velocities of FOF clusters in different mass 
intervals and found that the rms cluster velocity decreases 
with mass. Suhhonenko & Gramann (2003) studied the ve- 
locities of FOF clusters for different masses and radii. They 
showed that the effect of the cluster radius on v rms is similar 
to the effect of the cluster mass on v rms . The rms velocity of 
small FOF clusters is higher than the rms velocity of large 
massive clusters. 

We also plot the linear theory rms peculiar velocity 
of peaks, a p , defined by the linear radius R = Rl (dot- 
ted lines in Fig. 1). In this approximation, the rms pecu- 
liar velocity decreases as the cluster mass increases. For the 
mass M = lO 14 ft -1 M , the linear radius R L = 6.59ft -1 
Mpc and the rms velocity ct v {Rl) = 427 kms -1 . For 
M = 5 x 10 14 ft -1 M Q , we find that R L = 11.3ft -1 Mpc 
and <j p (Rl) = 382 kms -1 . These values for the rms veloci- 
ties are significantly lower than the rms velocities of clusters 



found from the simulation. This result was first obtained by 
Colberg et al. (2000), who compared the cluster velocities 
with the linear rms velocities of peaks at the radius R = Rl- 



4 VELOCITY DISTRIBUTION FOR THE 
CLUSTERS 

Fig. 2 and Fig.3 show the velocity distribution function for 
the DMAX clusters. We investigated the distribution of one- 
dimensional cluster peculiar velocities, (v x ,v y ,v z ), for differ- 
ent cluster masses. Fig. 2 shows the results for the radius 
R s = 1.5ft -1 Mpc and Fig.3 for the radius R s = 1.0ft -1 
Mpc. The probability density was estimated as the normal- 
ized number of clusters in the range vi ± Avi , as a function 
of vi. We used the value Avi = 100 kms -1 . We also show 
the Poisson error bars for the velocity distribution. 

We examined the velocity distribution for the clusters 
in three different mass intervals: (1) M = (5 x 10 12 — 

10 13 ) ft -1 M o (Fig. 2a and Fig. 3a), (2) M = (10 13 - 5 x 

10 14 ) ft -1 M o (Fig. 2b and Fig. 3b) and (3) M > 10 13 ft -1 M Q 
(Fig. 2c and Fig. 3c). The number of clusters, N c i, in the 
mass intervals (1) and (2) is shown in the Table 1. The num- 
ber of clusters with masses M > 10 13 ft -1 M Q is N cl = 1973 
and N c i = 1526 for the radius R s — 1.5ft -1 Mpc and 
R 3 = 1.0ft Mpc, respectively. We see that the velocity 
distribution function for different mass intervals is similar. 

For comparison, we plot in Fig. 2 the Gaussian distri- 
bution with the dispersion a\ given by o\ — a' p /\/3 = 
260 kms -1 . This dispersion is predicted by the linear the- 
ory for the scale R = 1.5ft -1 Mpc, after accounting for the 
finite size of the simulation box (see eq. 13). At the radius 
R = 1.5ft -1 Mpc, the rms velocity of peaks a p = 500 kms -1 
and, therefore, the parameter a' p = 450 kms -1 . Fig. 2 
demonstrates that the one- dimensional velocity distribution 
function for R s = 1.5ft -1 Mpc clusters is well approximated 
by a Gaussian. 

In Fig.3 we compare the velocity distribution of R s = 
1.0ft -1 Mpc clusters with the Gaussian distribution for 
(7i = 265 kms -1 . At the radius R = 1.0ft -1 Mpc, the 
rms velocity of peaks a p = 510 kms -1 and, therefore, 
a' = 460 kms -1 . We see that the velocity distribution of 
R s = 1.0ft -1 Mpc clusters is also similar to the Gaussian 
distribution. However, there are small deviations from the 
Gaussian distribution, especially for the clusters with masses 
M = (10 13 - 5 x 10 13 )ft -1 M Q . 

During the evolution larger and larger scales become 
non-linear and deviations from the Gaussian distribution de- 
velop. As demonstrated by Bachall, Gramann & Cen (1994) 
and Sheth and Diaferio (2001), evolution of the cluster pe- 
culiar velocities depends on the large-scale density field. The 
highest velocity clusters frequently originate in dense super- 
clusters. Fig. 3 demonstrates how the deviations from the 
Gaussian distribution start to develop. 

Fig. 4 shows the velocity distribution function for the 
FOF clusters determined by b = 0.2. We studied also the 
velocity distribution for the clusters defined by b — 0.15 and 
found similar results. The velocity distribution was deter- 
mined in different mass intervals. For comparison, we plot 
the Gaussian distribution function with the dispersion a 2 , 
as determined for the v x velocities in the simulation. The 
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Figure 2. Probability densities of one-dimensional cluster veloc- 
ity vi for limited mass ranges (denoted in each panel). Clusters 
are defined by the DMAX method with the radius R a = 1.5h~ 1 
Mpc. The solid, dot-dashed and dotted lines show the distribu- 
tion of the three Cartesian components (v x ,v y ,v z ), respectively. 
For clarity, the error bars are shown only for the v x distribu- 
tion. The heavy dashed line in each panel describes the Gaussian 
distribution function predicted by the linear theory for the scale 
R = 1.5ft" 1 Mpc (ai = 260 kms' 1 ). 
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Figure 3. As for the previous figure, but for the DMAX clusters 
defined with the radius R s = l.0h~ 1 Mpc. The heavy dashed 
line in each panel describes the Gaussian distribution function 
predicted by the linear theory for the scale R = 1.0h~ Mpc 
(<ri = 265 kms- 1 ). 



parameter a x — 266, 264 and 246 kms' 1 in the Fig. 4a, 
Fig. 4b and Fig. 4c, respectively. 

The velocity distribution of FOF clusters in different 
mass intervals is different. In Fig 4a and Fig. 4b, we see that 
the peculiar velocity distributions found for the clusters have 
stronger peculiar velocity tails than predicted by a Gaussian 
distribution. Deviations from the Gaussian distribution are 



larger for small low-mass clusters. This result for the FOF 
clusters was obtained also by Hamana et al. (2003), who 
studied the velocity distribution for clusters with masses 
M < 3.2x lO 13 ft _1 M0 (see Fig.5 in their study). We also in- 
vestigated the velocity distribution for massive clusters with 
M > 5 x 10 13 /i _1 M o (Fig. 4c). The velocity distribution of 
massive FOF clusters is similar to the Gaussian distribution. 
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Figure 4. As for the figures Fig. 2 and Fig. 3, but for the FOF clus- 
ters determined by b = 0.2. The heavy dashed line in each panel 
describes the Gaussian distribution function with the dispersion 
<r^ as computed for the v x velocities in the N-body simulation. 



5 EVOLUTION OF CLUSTER VELOCITIES 

We investigated the evolution of cluster velocities by using 
the Virgo simulation for the ACDM model at redshift z = 
10. 

We selected 509 clusters with masses M > 1.2 x 
10 14 /i _1 Mq at z = 0. The mean intercluster separation 
in this sample was d ci = 30ft -1 Mpc. The rms veloc- 
ity of the clusters was v rms — 494 kms - . The clusters 
were determined with the DM AX method, using the radius 
R s — 1.5ft. -1 Mpc. For each cluster, we found all particles 
within a 1.5ft -1 Mpc sphere around the centre of the cluster 



and determined the positions of these particles at z = 10. 
We define the initial centre of the cluster at z = 10 to be 
the barycentre of these particles. 

Are these initial centres of the clusters associated with 
peaks in the density field at z = 10? We studied the den- 
sity field in the Virgo simulation at z — 10. To identify the 
peaks in the density field, we used the DMAX method with 
the radius R s = 1.5ft -1 Mpc (comoving). The rms density 
contrast on the 256 3 grid was <ro = 0.31. For each peak, we 
determined the dimensionless height of a peak, v, and the 
peculiar velocity of a peak, v p . The parameter v was defined 
as 



1 

(To 



N 



(16) 



where Nd was the number of particles in the sphere of ra- 
dius R s around the centre of the peak and N was the mean 
number of particles in this sphere (see eq. 8) . The peculiar 
velocity of a peak was defined to be the mean velocity of par- 
ticles in the sphere of radius R a . The velocities of the peaks 
were scaled up to the values expected at z — according to 
the linear theory. In the linear regime, the peculiar velocity 
v ~ aDHf, where a is the scale factor and D is the linear 
growth factor. By taking into account the evolution of these 
four functions, we found that v(z = 0) = 2.42 v(z = 10). In 
the following we consider only peaks with v > 3. The rms 
peculiar velocity of these peaks was 516 kms - (scaled up 
to z = 0). 

We investigated the peaks within a 3.0ft _1 Mpc sphere 
around the initial centre of the clusters. We found that 413 
clusters (81 per cent) are associated with a peak or peaks in 
the density field at z = 10. Most of the remaining clusters 
can be associated with a peak of either a slightly lower height 
or at a slightly greater separation. For 205 clusters, we found 
one peak at the initial centre of the cluster. 144 clusters were 
associated with two peaks. For 64 clusters, we found three 
or more peaks. 

In Fig. 5 we compare the peculiar velocities of peaks at 
the initial centres of clusters with the final peculiar velocities 
of clusters at z = 0. In Fig. 5a we show the velocities for 
the clusters associated with one peak. Fig. 5b describes the 
results for the clusters associated with two peaks. We see 
the correlation between the peculiar velocities of clusters and 
their associated peaks. In Fig. 5a, the rms difference in (3-D) 
peculiar velocity between a cluster and its associated peak 
is 246 kms -1 . In Fig. 5b, the rms difference is 262 kms -1 . 

At z = 0, the density field is highly non-linear at the 
scale R = 1.5ft -1 Mpc. The rms density contrast is 5.03. 
We cannot consider the density maxima in this field as sep- 
arate objects, which move with linear speeds. During the 
evolution many density maxima merge and lose their iden- 
tity. Some clusters obtain higher velocities than the values 
expected according to the linear scaling, other clusters move 
at lower velocities. However, these effects compensate each 
other. Our results show that the rms peculiar velocity of den- 
sity maxima at z = is well approximated by the linear rms 
velocity expected on a given scale. The rms peculiar velocity 
of density maxima is similar to the rms peculiar velocity for 
all field points at the same scale. In other words, the rms 
peculiar velocity does not depend much on the location of 
spheres used to determine the velocities. In comparison with 
the density field, the velocity field is more heavily weighted 
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Figure 5. Peculiar velocities of peaks at the initial centre of the 
clusters, v p , as compared with the final peculiar velocities of clus- 
ters, v c i , at z = 0. The velocities of peaks are scaled up to values 
expected at z = according to the linear theory, (a) Velocities 
for the 205 clusters associated with one peak, (b) Velocities for 
the 144 clusters associated with two peaks. Different symbols in 
the panel (b) show the velocities for the different peaks. 



by modes with low values of the wavenumber k, i.e. large 
scales which are in the linear regime. 



6 SUMMARY AND DISCUSSION 

In this paper we have examined peculiar velocities of clusters 
predicted in the ACDM model. We analyzed the clusters in 
the Virgo simulation for the ACDM model with flo = 0.3, 
h — 0.7 and as = 0.9. To identify clusters in the simu- 
lation we used two methods: (1) the standard friends-of- 
friends (FOF) method and (2) the method, where the clus- 
ters are defined as the maxima of the smoothed density field 
(DMAX). We used a top-hat window with smoothing radii 
Rs = 1.5k' 1 Mpc and R s = 1.0ft -1 Mpc. The peculiar ve- 
locity of DMAX clusters was defined to be the mean peculiar 
velocity within a sphere of radius R s . 

We studied the rms peculiar velocity of clusters for dif- 
ferent cluster masses. We found that the relation between 
the rms peculiar velocity of a cluster and the mass of the 
cluster depends on the method that is used to determine 
the cluster masses and velocities. The rms peculiar veloc- 
ity of FOF clusters decreases as the mass of the clusters 
increases. The rms peculiar velocity of DMAX clusters is 



almost independent of the cluster mass and is well approxi- 
mated by the linear rms velocity of peaks of the density field 
smoothed at the radius R — Rs. We also studied the rms pe- 
culiar velocities of peaks, a p , for the linear smoothing radius 
R — Rl(M). In this approximation, the rms velocities are 
significantly lower than the rms velocities of clusters found 
in the simulation. 

We investigated the distribution functions of the cluster 
peculiar velocities for different cluster masses. The peculiar 
velocity distributions found for the low-mass FOF clusters 
have stronger tails than predicted by a Gaussian distribu- 
tion. But the velocity distribution of massive FOF clusters is 
similar to the Gaussian distribution. The velocity distribu- 
tion of DMAX clusters in different mass intervals is similar. 
The one-dimensional velocity distribution of R s = 1.5/i -1 
Mpc clusters is well approximated by a Gaussian. In the ve- 
locity distribution of R s = 1.0ft -1 Mpc clusters, we found 
small deviations from the Gaussian distribution. 

We also investigated the evolution of cluster velocities 
on the scale R s = 1.5/i _1 Mpc. We found a correlation be- 
tween the peculiar velocities of massive clusters at z = and 
the peculiar velocities of peaks at the initial centres of the 
clusters at z = 10. Density maxima at the scale R s = 1.5ft _1 
Mpc are not isolated objects that move with linear speeds. 

In this paper we introduced the DMAX method. This 
method operates with the smoothed density and velocity 
fields on a given scale Rs. We expect that the properties 
of DMAX clusters (e.g. their mass distribution, correlation 
function) can be determined analytically from the properties 
of initial density and velocity fields, at least, for the quasi- 
linear scales. The density and velocity distribution function 
for the maxima in the Gaussian density field was derived by 
Bardeen et al. (1986). Further study is needed to consider 
different properties of the DMAX clusters. 

We found that the rms velocities of DMAX clusters on 
a given scale R s are well described by the linear theory 
at R = R s (equation 4). Suhhonenko & Gramann (2003) 
studied the rms peculiar velocities of clusters in the rCDM 
model. They found that also in this model, the rms peculiar 
velocities of DMAX clusters are close to the linear theory ex- 
pectations. We can use the linear approximation to estimate 
rms peculiar velocities of clusters in different cosmological 
models. 

Fig. 6 shows the rms peculiar velocities of clusters esti- 
mated for different ACDM models using the approximation 
(4). Fig 6a shows the results for the radius R — 1.5ft _1 
Mpc and Fig. 6b for the radius R = 1.0ft -1 Mpc, respec- 
tively. We estimated the rms peculiar velocity of clusters 
in flat ACDM models, where the density parameter Q. m — 
0.2—0.35, the normalized Hubble constant ft = 0.65—0.7 and 
(jg = 0.7 — 0.9. The power spectrum of density fluctuations 
was chosen as P(k) — AkT 2 (k), where the transfer function 
T(k) was calculated using the fast Boltzmann code CMB- 
FAST developed by Seljak & Zaldarriaga (1996). The nor- 
malization constant, A, was chosen by fixing the value of as- 
On large scales, the power spectrum calculated by using the 
code CMBFAST is higher than the power spectrum given by 
equation (5) (for a fixed as; see, e.g., Fig. 2 by Gramann & 
Suhhonenko 2002). To determine the dimensionless growth 
rate /, we used the equations given by Peebles (1984). 

Consider the peculiar velocities of clusters for the 
smoothing radius R = 1.5ft _1 Mpc. If ft — 0.7 and as = 0.7, 
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Figure 6. The rms peculiar velocities of clusters predicted in 
different flat ACDM models, (a) The peculiar velocities are 
smoothed on the scale R = 1.5/1" 1 Mpc. (b) The smoothing ra- 
dius R = l.0h~ 1 Mpc. The solid lines show the results for h = 0.7 
and the dot-dashed lines for h = 0.65. The upper curves describe 
the velocities for as = 0.9 and lower curves for erg = 0.7. Dotted 
lines in panel (a) describe the rms peculiar velocity obtained by 
Giovanelli et al. (1998) for 24 observed clusters. They found that 
the rms cluster peculiar velocity is v rm s = 520 ± 139 kms -1 . 



the rms peculiar velocity of clusters is 403 kms -1 for 
fim = 0.2 and 426 kms -1 for Q, m = 0.35. For a higher 
density parameter Q. m , the function / is larger, but the am- 
plitude of the large-scale density fluctuations is smaller. If 
h = 0.7 and cr 8 = 0.9, the rms peculiar velocity of clusters 
is 518 kms -1 and 548 kms -1 for fl m — 0.2 and Q m — 0.35, 
respectively. 

For comparison, we show in Fig. 6a the rms peculiar ve- 
locity obtained by Giovanelli et al. (1998) for 24 observed 
clusters in the distance range ~ 10 and 90/i -1 Mpc. This 
comparison is very preliminary. The peculiar velocities for 
the observed galaxy clusters and for the simulated DMAX 
clusters are determined in different ways. Giovanelli et al. 
(1997) analysed the I-band Tully-Fisher (T-F) measure- 
ments for 782 spiral galaxies in the fields of 24 clusters (the 
SCI sample). Most of the galaxies used in the T-F analysis 
were within the Abell radius (Ra = 1.5fr -1 Mpc) around 
the cluster centre. The peculiar velocities for the clusters 
in this sample were studied by Giovanelli et al. (1998). In- 
dividual cluster T-F relations were referred to the average 
template relation to compute cluster peculiar velocities. Gio- 
vanelli et al. (1998) found that the rms one-dimensional clus- 



ter peculiar velocity in the SCI sample is 300 ± 80 kms -1 , 
which corresponds to the three-dimensional rms velocity 
Vrms = 520 ± 139 kms -1 . This number is in good agree- 
ment with previous estimates by Bahcall & Oh (1996) and 
by Watkins (1997), based on the SCI sample. This estimate 
is also in agreement with that determined for the SCII sam- 
ple (Dale et al. 1999b). The SCII sample is based on T-F 
measurements for 522 late-type galaxies in the fields of 52 
Abell clusters in the distance range ~ 50 to 200fr -1 Mpc. 
The distribution of these galaxies in the 35 Abell clusters 
was presented by Dale et al. (1999a). Dale et al. (1999b) 
studied the cluster peculiar velocities in the SCII sample 
and found that the rms one- dimensional cluster peculiar ve- 
locity in this sample is 341 ±93 kms -1 . This corresponds to 
the three-dimensional rms velocity 591 ± 161 kms -1 . 

Yoshikawa, Jing & Borner (2003) analyzed the velocity 
dispersion of galaxies, a ga i, and of dark matter, a dm, in- 
side clusters using hydrodynamical simulations. They found 
that the velocity dispersions of galaxies in less massive 
(~ 1O 13 M0) clusters are systematically lower than those of 
dark matter particles inside the clusters. We do not exactly 
know how the mean peculiar velocity of galaxies is related 
to the mean velocity of matter in the cluster. It is usually 
assumed that these velocities are similar. Further work is 
needed to study the peculiar velocities of galaxy clusters in 
the hydrodynamical simulations. 
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